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Abstract 

We examine the characteristics of the 7 r /2 pulse for simultaneously rotating 
two spin species of different gyromagnetic ratios with the same sign. For a 
7r/2 pulse using a rotating magnetic held, we derive the equation relating 
the frequency and strength of the pulse to the gyromagnetic ratios of the two 
particles and the strength of the constant holding held. For a 7 r /2 pulse using 
a linear oscillatory magnetic held, we obtain the solutions numerically, and 
compare them with the solutions for the rotating 7 r /2 pulse. Application of 
this analysis to the specihc case of rotating neutrons and ^He atoms simul¬ 
taneously with a 7r/2 pulse, proposed for a neutron electric dipole moment 
experiment, is also presented. 
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The 7 r /2 rotation is a commonly used technique in nuclear magnetic res¬ 
onance. A spin S with a gyromagnetic ratio 7 , pointing along a constant 
magnetic held Bq in the z direction, can be rotated into the x — y plane by 
applying a radio-frequency (rf) magnetic held B(t), called a 7 r /2 pulse, in the 
direction orthogonal to Bq. For a 7 r /2 pulse rotating in the x — y plane, B{t) 
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can be written as 


B^(t) = Bi sin {oJrft), By(t) = Bi cos {tJrft), (1) 

which represents a magnetic field with freqnency ojrf and an amplitnde Bi. 
While Urf is often chosen to be the same as the resonant Larmor freqnency, 
ojq = 7-00, off-resonance freqnencies conld also be nsed. As discnssed below, 
this capability to perform a 7r/2 rotation with off-resonance freqnencies makes 
it possible to rotate simnltaneonsly two spin species of different gyromagnetic 
ratios with a single 7r/2 pnlse of snitable freqnency and dnration. 

The need to rotate two different spin ^ecies with a single 7r/2 pnlse is 
relevant for the proposed experiment [H, y, 1^ at the Oak Ridge National 
Laboratory to search for the nentron electric dipole moment (EDM) nsing 
polarized nltra-cold nentrons stored in a snperfinid helinm cell containing 
polarized ^He acting as a co-magnetometer and a spin-analyzer. A non-zero 
nentron EDM can be observed by measnring the difference of the precession 
freqnencies of nentrons when electric and magnetic fields are aligned or anti¬ 
aligned. The initial polarization directions for nentron and ^He spins are 
parallel to a constant Bq field pointing along the 2 axis. To measnre the 
precession freqnency of nentrons relative to that of ^He, a single 7r/2 pnlse 
will be applied to rotate both spins into the x — y plane. A snpercondncting 
qnantnm interference devices (SQUID) will be nsed to measnre the precession 
freqnency of ^He j^, and the relative precession between nentrons and ^He 
can be determined by the rate of the absorption reaction, n -|- ^He —t p -|- ^H, 
which depends sensitively on the relative spin orientations of nentrons and 
^He [^. From the measnrements of the precession of ^He and the absorption 
reaction of n-^He, the precession freqnency of nentrons can be determined. 
An alternative method, the dressed spin techniqne, conld also be nsed to 
measnre the nentron EDM by applying an additional dressing field [H, [^. 

It is not evident that two spin species with different gyromagnetic ratios 
conld be rotated simnltaneonsly by a 7r/2 pnlse. As shown in this paper, this 
is only possible if the ratio of the magnetic moments of the two spin species 
falls within a certain range. Vasserman et al. ntilized a method to simnl¬ 
taneonsly rotate the spins of electron and positron by 90° to measnre their 
anomalons magnetic moments for a test of CRT invariance . This method 
is only applicable for particles with identical magnitnde of gyromagnetic ra¬ 
tios. Recently, de Lange et al. pnblished a new method to manipulate two 
spin species using the spin echo technique [sj. However, this method requires 
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a spin bath which does not exist in many experiments such as the neutron 
EDM experiment. It is necessary to consider other methods which do not 
require a spin bath. The purpose of this paper is to discuss a method which 
is generally applicable to experiments requiring a simultaneous 7r/2 rotation 
for two different spin species. 

The rest of this paper is organized as follows. First, we derive the equa¬ 
tions for calculating Urj and r as a function of Bq, Bi, and the gyromagnetic 
ratios of two different spin species. Second, we explore the range of these 
parameters and the relation between them. Third, we present the numerical 
solutions and comment on the characteristics of the solutions. Finally, we 
also consider the case for linear rf magnetic helds, which are more readily 
implemented than rotating helds, and compare the solutions with those of 
rotating rf helds. The application of this study to the specihc case of the 
neutron EDM experiment involving neutron and ^He will also be presented. 
Some initial results obtained by one of the coauthors (JC^ for the case of 
rotating rf helds were presented in an unpublished report [2(]. 

In general, the spin motion of polarized particles moving in an exter¬ 
nal electromagnetic held can be described by the Bargmann-Michel-Telegdi 
equation [^. However, for non-relativistic particles, the dynamics of a spin 
in a magnetic held system can be described by the Bloch equation (lol |: 

^ = S X (7B(()). (2) 

Here we consider a system of polarized particles moving slowly in a cell with a 
static magnetic held, Bq, parallel to the spin orientation of polarized particles 
along the z-axis. Upon the application of the rotating held given in Eq. [H the 
ehect on the spin direction of the particles can be conveniently described by 
transforming to a frame rotating at Urf along 5. In this rotating frame with 
axes dehned by the x',y', and z' as shown in Fig. [H the held B{t) becomes 
static with a magnitude Bi pointing along the y' axis. The static Bq held in 
the z' axis is changed to Bq — in this rotating frame and the spins of 

the particles will precess about the magnetic held Btot given as 

Btot = [Bq - ^)z' + Biy' (3) 
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with an ehective Larmor frequency 

= 7^(50 - + BI (4) 
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Figure 1: (Color figure) Schematic plot of spin rotation using a rotating rf field Bi{t) 
in the x — y plane with frequency uJrf- The holding field Bq is along z-axis. In a frame 
rotating at uJrf along z, the effective field is Btot = {Bq — + Biy'. S'(O) is the initial 

spin orientation and S{t) is the spin orientation after applying a 7r/2 pulse, tp is the angle 
between Btot and y' and 4> the angle between S(t) and x'. 


The angle -0 between the directions of Btot and y' satishes 


(5) 


It is convenient to dehne another coordinate system {x' , y', z') obtained 
from a rotation along the x' axis by an angle 'ip snch that y is along the 
direction of Btot'- 


- ^ - 

X 


■ 1 0 o' 


x' 

^ II 

y 

= 

0 cos ip sin pj 


y' 

^11 

Z 


0 — sin -0 cos ip 


z' 


The spin of the particle is initially along the i'-axis. At a later time t, it can 
be expressed as 


S{t) = S^pt) + S^{t), 


(7) 
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where S\\(t) is the component parallel to Btot and is independent of time, 

^ll(t) = l^lsinV'/. (8) 

S± is the component perpendicular to Btot and will rotate about Btot so that 

S±{t) = {\S\ cos'll)) {cos {uQt)z —sm{uQt)x ). ( 9 ) 

Using Eqs. E]- [HI S{t) can be expressed in terms of the coordinates in the 
rotating frame as 

S{t) = S\\{t) + S^{t) 

= 1^1 (sin+ ( 1^1 cos V’)(cos (ci;Qf)z — sin (cjQt)^ ) 

= I*?! [— cos-^ sin(a;Qf)x' + sin V' cos'0(1 — cosoj^t)!)' 

+ (sin^ 0 + cos^ 0 cos . ( 10 ) 

For an rf pulse to rotate the spin from the vertical axis to the horizontal 
plane after a duration r, the component of S{t) along z' must vanish: 

S{t) ■ z' = sin^ 0 + cos^ 0 cos oj'qT = 0 . ( 11 ) 

We can readily obtain the solutions for and r to rotate two spin 
species of different gyromagnetic ratios by 7 r /2 simultaneously. Due to the 
property of the cos function, the sign of Uq does not change the solution of 
r; this implies the solution of r is suitable for both positive and negative 
Uq as well as 7. Denoting the gyromagnetic ratios for the two spin species 
as 7 i and 72 , and their Larmor frequencies as cuop = IiBq and Uo ,2 = 72 -Bo; 
the requirement is to have the same duration, r, with a given frequency to 
simultaneously rotate both species by 7 r/ 2 . Therefore, we obtain: 

cos“ 0 -tan 2 ( 0 i)) cos"^(-tan^( 02 )) 

^ ~ l77 i ~ l77 i 

1 ^ 0 , 11 1 ^ 0 , 2 ! 

where oJq ^ and 0* are the effective Larmor frequency and 0 angle for species 
i in Eq. [Hand Eq. [5l respectively. The range of cos“^(—tan^(0)) is from 0 
to TT in order to keep r positive. 

The solutions for ujrf must also satisfy the following constraint: 


— 1 < tan0j < 1, 
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i.e. 


implying 


-Dl 


^rf 

liBo 


< 1 , 


Bl ^ OJrf 

Bq 7iBo 




(13) 


(14) 


It is also noted that after the application of a 7 r /2 pulse, the spin directions 
of the two species on the x — y plane would in general be different. Using 
Eqs. [To] and mi the azimuthal angle 0* of the spin for species i after Q.'nj‘2 
rotation is given as 


= tan 


sin 0 i cos 0,(1 — coscug jT) , 
— cos ipi sin Uq {t 


= — tan 


-1/ 


ycotU/j) 


(15) 


Therefore, the spins of the two species are usually not aligned (0i ^ 02 ) right 
after the 7r/2 pulse. 

Equation [12] can be solved numerically to hnd Urf for given values of 
Bq, Bi, and y*. Because the gyromagnetic ratios of neutron and ^He are both 
negative, we focus on the case of both spin species having gyromagnetic ratios 
of the same sign in this paper. We note that Eq. [T2] is also applicable for the 
case of two spin species having gyromagnetic ratios of opposite signs. For 
convenience, we define dimensionless parameters if both gyromagnetic ratios 
are in the same sign: R = ^, x = ^ and y = Fig- E^a) shows the 

solutions for y as a function of x for three different values of R. We only 
consider the case for i? > 1, since the R < 1 case simply corresponds to 
interchanging the two spin species. We first note that for a given value of 
R [R = 71 / 72 ), a solution for Urf can be found only within a certain range of 
X. This range becomes narrower as R increases, and above certain value of 
R, there is no longer a solution for Urf- The domain in R versus x for which 
a solution for Urf exists is shown as the red region in Fig. E^b). 

Figure E]^a) also shows that the solutions for y at different values of R 
have very similar shapes. It is interesting to study the phenomena when R 
approaches 1 (but not equal to 1). Inserting R ^ 1 + 6R into Eq. [12] and 
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ignoring higher-order terms in hi?, the solntion for y versus x when i? —)■ 1 is 

cos“^(a;^(—1 -I- yY) 

2{x‘^y{l-y)) /l + a;2(l -?/)2 
"" H-a;2(l-2/)y l-a;2(l-|/)2’ ^ 

shown as the black-dashed curve in Fig. [2](a). However, when i? = 1, corre¬ 
sponding to the degenerate case of a 7r/2 rotation for only a single species, 
the solution is given by Eq. [141 shown as the yellow region in Fig. [2]^a). This 
abrupt change from a curve to a band for the solution of Eq. [12] reflects the 
degeneracy occurring at i? = 1. 

The solution for the specihc case of rotating neutron and ^He simultane¬ 
ously with a 7r/2 pulse, relevant for a proposed neutron EDM experiment [^, 
is shown in Fig. [3l The black-dashed curve in Fig. [3](a) gives the solution for 
y versus x aX R = 'js/'jn = 1-1121, where n and 3 represent neutron and ^He, 
respectively. Figure [31^b) and (c) also show solutions of the duration r from 
Eq. [T2]and the angular difference of spins of two species after a 7r/2 rotation, 
~ 03, from Eq. [TS] Also shown in Fig. [3] are the solutions for 
linear rf fields, to be discussed next. 

The study so far assumes a rotating rf held. However, a linear rf held is 
commonly utilized in experiments due to its simplicity to implement. There¬ 
fore, we extend our study to the case of a linear rf held, which can be de¬ 
composed in terms of two rotating components: 

Brf{t) = 2Bicos{urft)y (17) 

= Bi{sm{Urft)x + COs{Ur ft)y) + Bi{— sm{urft)x + COs{Urft)y) 

The hrst component corresponds to the rotating held we have considered, 
which is stationary in the rotating frame. The second component rotates 
with a frequency of —2urf in the rotating frame. To assess the ehect of 
this high frequency term, we have solved the time dependence of S numeri¬ 
cally using the Bloch equation (Eq. [2]). We have considered the specihc case 
for neutron and ^He with the gyromagnetic ratio —18.32472 and —20.37895 
Hz/mG, respectively, and with Bq at 10 mG. The Runge-Kutta method is 
applied and the time step of the simulation is At = 10“® sec. In order to 
determine a single 7i/2 pulse for both neutron and ^He, the following algo¬ 
rithm is applied. Initially, both neutron and ^He spins are along the i-axis. 
Using the numerical simulation of the Bloch equation for a given x = ^ and 
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y — —derive the duration r when the neutron spin is rotated into 
the X — y plane, i.e., the ^-component of the neutron spin becomes zero. In 
practical terms we derive the duration when the z-component of the neutron 
spin changes its sign, from a positive value to a negative value. Although 
the spin is wobbling, we just use the hrst point when the neutron spin is 
rotated into the x — y plane. Then, at r, we consider the difference of the 
i-component between the spins of neutron and ^He. For a given x = we 

vary y = —— until the difference of the z-component between the spins 
of neutron and ^He changes sign. Using this algorithm, we can derive the 
solutions of the linear rf helds for the simultaneous 7r/2 rotation of neutron 
and ^He. 

Fig. m shows an example of a 7 r /2 pulse for neutron and ^He using a linear 
rf held. The vertical components of the spins of the two species as a function 
of time are shown as red-solid and green-dotted curves for neutron and ^He, 
respectively. The oscillatory pattern in these two curves has a frequency of 
~ 2urf. Fig. m shows that the both spins can be simultaneously rotated to 
the X — y plane with the linear rf held. 

In Fig. [3]^a), the solutions using the Bloch equation are shown for both 
rotating and linear rf helds. Signihcant diherence between the two cases is 
observed for smaller values of BqIBi. At large values of Bq/Bi, the ehect 
of the high-frequency counter rotating term could be neglected. The time 
duration r and the angular diherence after 7r/2 pulses for i? = 73 / 7 ^ are also 
shown in Fig. [3t^b) and (c). They show that the results of the time duration 
in Eq. [12] and the angular diherence in Eq. [15] are consistent with those of 
numerical simulation using Bloch equation for rotating rf helds. The results 
for linear rf helds follow those for rotating rf helds, but with an oscillatory 
pattern superimposed. 

In summary, we have studied the solutions for 7 r /2 pulses for two spin 
species of the same-sign gyromagnetic ratios using rotating and linear rf 
helds. The characteristics of the solutions are presented. For a specihc 
experiment, the selection of the optimal values of Bo, Bi could depend on 
various considerations. This study provides the solution for Urf, once the 
values of Bo,Bi are chosen. The application of this work to a neutron EDM 
experiment [2| is also discussed. Finally, this study can be extended to the 
case of two spin species having gyromagnetic ratios of opposite signs. 

We gratefully acknowledge valuable discussions with Bradley W. Filip- 
pone and Riccardo Schmid. This work was supported by the U.S. National 




Science Foundation and the Department of Energy. 
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x=Bo/B^ 


Figure 2: (Color figure) Solutions of the 7 r /2 rotation for two spin species, (a) y = 
2a;r//(cLio,i + ^ 0 , 2 ) versus x = Bq/Bi for different R = 71 / 72 - The yellow region shows the 
constraint of the 7 r /2 rotation for a single spin species according to Eq. [Ml (b) i? = 71/72 
versus x = Bq/Bi. All values in the red region can be used for 7 r /2 pulses. 
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x=Bq/Bi 


Figure 3: (Color figure) Solutions of the 7r/2 pulses for neutron and ^He, i.e., R = 

The holding field Bq is 10 mG along the z-axis. (a) y = 2uJrf/{ ujq^i + 0 ^ 0 , 2 ) versus x = 
Bq/Bi. (b) Time duration r versus Bq/Bi. ( c ) The angular difference, between 

spins of neutron and ^He after the 7r/2 pulses. The black-dashed curves and the red-solid 
curves are solutions for rotating rf fields using the analytical equations (Eqs. [T^] and [T5]) 
and the Bloch equation, respectively. The results from these two approaches are identical. 
The green-dotted curves are obtained using Bloch equation simulation with linear rf fields. 
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Figure 4: (Color figure) The time dependence of Sz for neutron and ^He with a linear rf 
field. The red-solid and green-dotted curves are neutron and ^He, respectively. Bq is 10 
mG along the £-axis, and the black-dashed curve shows the amplitude and frequency of 
the linear rf field, Brf(t), defined in Eq. [T71 (x = Bq/Bi = 17, y = 2 “ 0.989596) 
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